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Abstract 

Dynamic holography has been demonstrated as a method for correcting aberrations in space 
deployable optics, and can also be used to achieve high-resolution beam steering in the same 
environment. In this paper, we consider some of the factors affecting the efficiency of these devices. 
Specifically, the effect on the efficiency of a highly collimated beam from the number of discrete phase 
steps per period is considered for a blazed thin film beam steering grating. The effect of the number of 
discrete phase steps per period on steering resolution is also considered. We also present some result of 
Finite-Difference Time-Domain (FDTD) calculations of light propagating through liquid crystal "blazed” 
gratings. Liquid crystal gratings are shown to spatially modulate both the phase and amplitude of the 
propagating light. 


1. Introduction 

As NASA develops a need to transmit greater amounts of data over long distances, there is a 
need to develop new technologies capable of such performance. 1 Current S- and X-band RF technology 
provides only limited room for increasing the data communication rate. ' The use of K a -band technology 
may increase the usable photon flux by a few orders of magnitude, but even that may not be enough to 
permit, for example, high-speed, long-distance communication of high -resolution image data. 4 In 
addition, increasing the RF transmitter power level does not get past the bandwidth limitations of the 
carrier frequency. For that reason, NASA is exploring the use of shorter wavelengths in the near IR, 
where the greater directivity of transmitting beams can increase the photon flux and bandwidth. 

However, this potential benefit is not without cost, for it involves trading problems with existing RF 
technology for a new set of problems. Among these new problems is the necessity for some means of 
high-resolution optical beam steering in order to maintain connection of the more directional signal beam 
with the distant receiver. Because of the anticipated demand for higher data rates at all transmission 
distances, special steps must be taken in order to ensure a high signal concentration at the receiving end. 
hi order to maintain the intensity of a propagating laser beam after long propagation distances, the beam 
must possess a divergence of a few microradians or less. Along with this directivity requirement comes 
the need for a steering capability of even greater accuracy, i.e., sub-microradian beam pointing 
technology. 5 
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The most mature beam steering devices achieve deflection through mechanical means. Although 
mechanical pointing technology is being considered for use in deep-space optical communications, 5 7 
such devices may not provide the needed precision. 

Liquid crystal devices have already been considered for use in beam steering applications, 8111 as well 
as for the closely related function of tip-tilt correction of wave fronts. 11 

In this paper we discuss factors influencing the range, resolution, and number of steering angles 
available from a "phase-only’' liquid crystal spatial light modulator (LC-SLM). Computational simula- 
tions are used to quantify the effect of phase stair-stepping on beam steering (i.e., diffractive) perform- 
ance. Simulation data is also presented to quantify the impact of beam truncation by the diffracting 
aperture and surface flatness. 

These devices can be designed at or between either of tw o extremes: large pixels producing stair- 
stepped phase modulation or small pixels whose director distortions produce smoothed phase profiles. 
Here, we discuss the large pixel variety. The effects of director distortion and finite cell thickness will 
be considered in later publications. 


2. Background 

The LC-SLM consists of two glass plates sandwiching a thin (~10 pm or less) layer of liquid 
crystalline material. One substrate is coated w ith a single common electrode, usually held at zero volts. 
The other glass substrate is patterned with a number of uniformly sized and spaced discrete electrodes. 

If suitably addressed, the electrodes of the LC-SLM selectively reorient the liquid crystal to form a 
one-dimensional phase grating. This grating is used to impose a periodic "blazed” phase profile onto 
the incident light (fig. 1). In beam steering or tilt-tip correction applications, the goal is to provide the 
equivalent of a linear phase ramp with its lit phase degeneracies removed (fig. 2). Operated in this 
fashion, the LC-SLM can be used as a means of redirecting the path of monochromatic light in one 
dimension. Two-dimensional beam steering or tip-tilt correction can be obtained by either cascading two 
one-dimensional LC-SLM’s, or by encoding onto a single device a two-dimensional array of individually 
addressable electrodes. 

The use of an LC-SLM in this manner has carried the label of Optical Phased Array (OPA) 
technology. 8 ’ 9 In one implementation of liquid cry stal OPA technology, 10 every rt electrode is linked 
together in order to reduce the number of off-device connections. This approach is taken primarily to 
make more practical the fabrication of large apertures containing a large number of electrodes. However, 
that approach also limits available steering angles to a discrete number of angularly separated deflections, 
hi order to completely cover an entire region of interest, cascading with a second fine-tuning stage is 
required. Technology does exist for connecting and individually addressing each of a moderately large 
number of (more than one thousand) electrodes. This latter form of LC-OPA technology is the subject of 
this paper. 

Prior analysis of the OPA, for which identical stair-stepped phase approximations of ideal blaze 
periods are repeated across the entire device, has yielded a theoretical maximum beam steering efficiency 
of: 8 


tj= sine 2 {tiI N) (1) 

where N is the number of stair-stepped phase segments in each identical blaze profile. It is assumed that 
all modulo-27r resets of the original ideal phase profile fall between electrodes. Below we will show that 
this result extends to the more general case of resets falling anywhere in the LC-SLM, even in the middle 
of electrodes, if the phase retardation produced by those electrodes is properly chosen. Hence, it is 
predicted that individually addressable lines will permit continuous steering from this stair-stepped LC- 
SLM device, without any penalty beyond that imposed by eq. (1). In other words. 
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this device does not, in theory, possess any pixel-related limit to angular resolution. This potential for 
achieving continuous steering from liquid crystal OPA technology has not, to our knowledge, been 
previously analyzed and discussed. 

Experimental measurement of the beam steering efficiency of one OPA-like liquid crystal 
grating was shown to fall short of that predicted by eq. (I). 11 The cause of that discrepancy is one of the 
motivations for this paper. 

As was mentioned, these devices can be implemented at or between two extremes. At one extreme 
the LC-SLM pixels or stripes (one for each electrode) consist of electrodes large enough that 
the phase profile exhibits stair-stepping. Fringe fields and liquid crystal elasticity at the edges of the 
electrodes (i.e., at the boundaries between adjacent pixels or stripes) cause some smoothing of this stair- 
stepping, but the area occupied by those behaviors is relatively insignificant relative to the much larger 
area of the electrodes as a whole. As a result, the device exhibits relatively sharp phase "resets,” closely 
approximating that of an "ideal” blazed phase grating (fig. 3a). At the other extreme, the pixels or stripes 
would be small enough that fringe fields and liquid crystal elasticity completely smooth out the phase 
profile, preventing the appearance of sharp phase "resets” (fig. 3b). This paper will be restricted to the 
former extreme. 


3. Methodology 

Assuming the thin-layer model, diffractive behavior can be obtained by means of a numerical 
evaluation of an appropriate scalar diffraction integral, such as the Integral Theorem of Helmholtz and 
Kirchhoff: 12 


y/{r) = j- J_[ s ,n • [y/{r')V'G - GV>(r')] dS' 


( 2 ) 


where ^is a complex scalar field representing light emerging from the grating layer, r' = (x'ji'^') is a 
location on the exit surface S” of the diffracting object, r = (xj/,z) is a location of some diffraction 
observation, n is the local normal to the surface of the diffracting object, and G is the appropriate Green’s 
function. For a one-dimensional grating such as is central to this discussion, the above equation is easily 
recast as a line integral: 1 ’ 


y/(r) = - 


cxp(-//r/4) 
(8 kx) m 


cxp( ikR ) 


R 


1/2 


d y/{x')+ l -^-y/{x') 

/VII 


\dx ' 


(3) 


where for a one-dimensional planar diffracting aperture R t = |r-r j, r = (x.y.O). r' = (x'.O.O). k is the 
magnitude of the wave vector, and the aperture extends over the range x<x\-d.d\. 

The two bracketed terms in eq. (3) are direct descendants of the two Rayleigh-Sommerfeld 
components of eq. (2). In general, if paraxial calculations are carried out for feature sizes within the 
aperture and propagation distance both much larger than the wavelength, the two Rayleigh-Sommerfeld 
components and the Helmholtz-Kirchhoff average of the two do not produce significantly different 
results, in particular for small deflection angles. The computations in this study evaluate the second term. 
For a Gaussian profile at the minimum beam waist, that term integrates the one-dimensional profile of the 
electric field of the layer-processed beam: 


y/{x') = Eq exp 


»\2 


-(*') 

w 0 


exp[/^(x')] 


(4) 
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where w 0 is the minimum beam waist radius and x' is the in-plane distance from the center of the beam, 
hi order to approximate thin layer phase gratings, the phase term, cxp|//7(.v')| can be tailored to the desired 
grating profile. For an unhindered Gaussian beam, (fix) = 0. 

Since the application of interest is optical communication to and from deep space probes, it is 
assumed that the propagation distance [/?,= |r-r' | in eq. (3)] is large enough to qualify the diffractive 
behavior as "far field” (beyond some propagation distance, the angular dependence of diffraction will be 
independent of that distance). 

All thin-layer approximated scalar diffraction calculations below were performed for a A = 1 pm 
Gaussian beam with a minimum waist of 25 cm coinciding with the "very thin” LC-SLM. All results are 
presented in the form of Strehl ratios. The Strehl ratio is defined here as the ratio of a particular diffracted 
intensity (field squared) to the peak intensity of an unhindered beam at the same propagation distance. 
That reference point was computed using the same Gaussian source, for an aperture diameter equal to 
8wVj. For calculations of the Strehl ratio, in which one intensity is divided by another, the value of E 0 is 
not important. 


4. Results 

The above calculation method was utilized to gauge the impact of beam truncation, surface flatness, 
and phase stair-stepping on the far field diffraction performance of large-pixel LC-SLM's. Since it is 
desirable to provide the maximum signal intensity to the receiver, it is important to know how each of 
these factors will influence the design of free-space optical communication components such as the beam- 
steering LC-SLM. 


A. Beam Truncation 

It is well known that if an aperture truncates a beam of light directed at it, the aperture not only blocks 
some of the incident light but also redirects some of the transmitted light into higher diffraction orders. 
This will also affect the Strehl ratio at the center of the propagating beam. For example, the simulated 
impact of beam truncation on a 25-cm-waist beam with X = 1 pm is shown in fig. 4. If the aperture's 
width is four times the minimum beam waist (which is assumed to occur in the aperture "plane”), the 
Strehl ratio measured at the center of the propagating beam in the far field is still a healthy 0.99. If the 
aperture size is reduced to three times the minimum waist, the Strehl ratio is reduced to 0.94, and if the 
aperture is further reduced to twice the minimum waist, to 0.72. Such losses would set an upper limit for 
performance, subject to other sources of loss. 


B. Imperfect Optical Surfaces 

If an optical surface in a beam steering device is less than perfect, those aberrations also reduce the 
amount of energy delivered to the intended direction. Optical surface imperfections are commonly 
decomposed into constituent aberration modes, which are then related to image-related effects. Liquid 
crystal devices employ "flat” manufactured transmissive and reflective glass substrates whose surfaces 
are perhaps better represented as a Fourier sum of surface undulation terms, each term of a different 
spatial frequency. One can then analyze the optical effect of each frequency component and determine 
which spatial frequency components, and how much of each, can be tolerated. The phase profile of such 
phase undulations, concave at the center of the aperture diameter, can be described by: 
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($(x) = 2 jiA cos 


2 m 


1 x x midpoint ^ 


V 


D 


j 


( 5 ) 


where A is one-half of the peak-to-valley amplitude of the undulation, x is the surface location, x nml|ioim is 
the center of the aperture, and D is the width of the aperture. This factor was substituted into eq. (4) to 
obtain the thin-layer approximation of the field transmitted through or reflected from the surface. For a 
reflective surface, the phase (j) in eq. (5) is multiplied by two, accounting for the total optical path offset. 

The effect of surface undulations can be computed for the case of light reflected from a nominally flat 
mirror. Assuming that the mirror size is four times the minimum waist (which coincides with the mirror 
surface), the effect of a surface undulation whose period equals the aperture size can be seen in 
fig. 5. Those results were obtained for the case of an undulation centered on the mirror, concave at the 
center. Even if the peak-to-valley undulation amplitude is a very fine 2/20. the resulting Strehl ratio at 
the center of the far field is reduced to 0.94. Performance for larger undulation amplitudes is noticeably 
worse (fig. 6). As the undulation frequency increases beyond the range of fig. 6, light will continue to 
be diverted into other diffraction orders until the undulations become smaller than a wavelength, at which 
point the surface will again appear to be flat. 

Of course, a glass or reflective silicon substrate employed in an LC-SLM could possess a range of 
undulation frequencies resulting from the manufacturing process, so no one spatial frequency will 
dominate. Nevertheless, it is clear that manufacturing imperfections will play a role and should be 
considered. 

This last point provides reason for concern. If space probes are to be equipped with meter-class 
optical apertures, it will not likely be feasible to launch glass apertures of thickness sufficient to provide 
stable, precise optical surfaces. Instead, optical components will likely be of lower mass, and thus not as 
rigid, and more prone to significant surface imperfections, hi order to make use of optical surfaces with 
substantial imperfections, it may be necessary to perform wavefront correction close to the source of the 
imperfections. It is not inconceivable that a two-dimensional LC-SLM could provide that function while 
simultaneously steering the beam. 


C. Phase Stair-stepping 

We now return to the issue of phase stair-stepping, which is the product of the segmented design of 
LC-SLM electrodes. As shown in figs. 2 and 3a, the principle behind OPA technology involves removal 
of 2/t degeneracies from a linear phase ramp prescribing the beam deflection angle, and then utilizing 
that blazed phase profile to assign phase values to the OPA's segments. This process is simple and 
straightforward if all of the "resets” fall between segments, as shown in fig. 3 a. Flowever, as was 
mentioned earlier, that approach permits only a finite number of discrete deflection angles, thereby 
placing a limitation on beam steering angular resolution. 

The angular resolution, as the term is used here, is not the same as the ability of a beam steering 
device to resolve spots. The maximum number of resolvable spots obtainable from a beam deflecting 
device has been defined as the device's maximum deflection angle divided by the beam's "diffraction 
limited angle.” 14 For a Gaussian beam that angle is defined as the divergence angle. If a beam steering 
device is engineered to provide deflection to only those separated angles, it will not be capable of 
continuously covering each point in its angular range with a signal of uniform intensity (see fig. 7). A 
device may be capable of resolving the maximum number of spots within its angular range, but unless 
it can steer to angles between the locations of those spots, its angular resolution is limited to the angle 
between those spots. In order to permit very high signal intensities everywhere within the device's 
angular range, angular steering resolution must be finer than spot resolution. 

Achieving sub-spot angular resolution (i.e., sub-microradian resolution, for deep space optical 
communications using meter-class transmitting apertures) may at times require that the original linear 
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phase ramp produce resets which do not fall between OPA segments (fig. 8). It is important that any such 
OPA segments produce the appropriate phase retardation level. If the phase of those segments is 
incorrectly assigned, diffraction performance (i.e., beam steering efficiency) will suffer. 

One obvious solution is to assign each of these "interperiod” segments a phase equal to the average 
value of the ideal blazed phase profile within the segment, as shown in fig. 9. This approach can then be 
described as applying two operations to the original linear phase ramp, namely modulo-2 ;r and area 
averaging, in that order. Those two operations can be carried out in reverse order, producing interperiod 
phase values as shown in fig. 10. 

The difference between these two approaches is illustrated in fig. 1 1 . For instance, applying the first 
approach to a linear phase ramp for which the modulo-2 ablaze width is equal to the width of 2 and 1/3 
electrode segments, the resulting phase profile will appear as shown in fig. 11a. The phase profile using 
the second approach appears in fig. 1 lb. It is only when the blaze width equals an integer number of 
electrode segments that the stair-stepped phase profile obtained using either approach will produce 
identical neighboring stair-stepped blazes as are found in fig. 3 a. 

Diffraction simulations were performed for grating phase profiles obtained using each approach. 

In each simulation, stair-step segments of 100 pm were used to define the stair-stepped phase profiles 
covering a lm wide one-dimensional aperture. Each simulation was carried out over a range of deflection 
angles, each of which is characterized by its own ideal linear phase ramp. Each linear phase ramp 
corresponds to a number of electrode segments per ideal blaze period, not necessarily an integer. At 
each value of electrodes per period, the diffracted intensity at the intended far field deflection angle was 
divided by that produced by the corresponding unbroken linear phase ramp to produce the Strehl ratio 
data appearing in fig. 12. 

For both methods of obtained stair-stepped phase profiles, the Strehl ratio at integer values of 
electrodes per period matches that obtained from eq. (1). For non -integer values, the first approach 
falls well below those values, while the second approach does not. Utilizing the second method, high- 
resolution steering and tip-tilt correction is possible from the LC-SLM. Performance does begin to 
falter at and below eight electrodes per period, as has been noted. s " Although it might appear that the 
"irregular periodicity” (see fig. 1 lb) of the more desirable second stair-step approach might violate the 
grating condition of its underlying blazed phase profile, such is not the case. The center of each of the 
stair-step segments does fall on the underlying linear phase ramp, and on its blazed-phase equivalent. 


D. Liquid Crystal Layer Simulations 

The diffraction simulations above employ approximations of the fields emerging from thin layer 
liquid crystal gratings. Those computations do not take into account optical anisotropy and spatial 
variations of the liquid crystal director required to produce diffraction gratings. In a real device, 
application of appropriate voltages to the closely spaced electrodes might lead to a director configuration 
which looks somewhat like fig. 1. In order to increase the level of accuracy of diffraction simulations, 
one must obtain more accurate surface fields than were utilized above. Direct analytical solution for the 
influence of such a structure on a propagating beam is at best very difficult. However, numerical methods 
such as the Finite-Difference Time-Domain (FDTD) are capable of providing a very accurate solution 
for such a problem. Such analysis may provide insight into the impact of anisotropy and inhomogeneity 
on diffractive performance. 

An illustrative two-dimensional simulation of those effects was conducted for this study. An 
idealized director configuration (fig. 13) was used as the repeated periodic unit of a hypothetical second- 
order liquid crystal grating designed to steer to 20 degrees. Segmentation was not encoded into the 
structure at this point, so as to promote the visibility of anisotropy and inhomogeneity -induced effects. 

The finite grating consisted of about 200 periods. The near field simulation employed FDTD to obtain 
the response of the idealized liquid crystal grating to an out-of-plane H polarized incident beam of 
Gaussian profile. 
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After obtaining a steady-state computed result, a snapshot was obtained of the fields in and beyond 
the grating. A section of this snapshot is shown in fig. 14, with the grating periods artificially outlined. 
While the grating appears to have imparted some deflection onto the transmitted fields, 
there is some disorganization in those fields, which may indicate reduced diffraction efficiency. The 
disorganization can be separated into amplitude and phase effects as shown in figs. 15 and 16. It appears 
from fig. 16 that the “phase-only” grating also modulates the transmitted intensity, which degrades the 
diffraction efficiency below the amount predicted by some imperfect phase modulation alone. 

A more complete study of FDTD simulations of segmented liquid crystal beam steering gratings will 
be carried out in the near future. 


5. Conclusions 

Implementation of optical communications with deep space probes will require very precise sub- 
microradian steering of the communication beam in order to maximize the received power. In addition, 
the received wavefronts may also require further correction, some component of which may be tip-tilt. 

We have analyzed the possibility of using liquid crystal spatial light modulators (LC-SLMs) to meet these 
needs. Three factors influencing the beam steering performance of these segmented liquid crystal devices 
have been analyzed: beam truncation, surface flatness, and phase stair-stepping. 

If transmitting a laser beam through an aperture, any truncation of the beam by the aperture detracts 
from the intensity arriving at a distant receiver. This effect begins to be noticed for apertures less than 
four times the waist of the beam passing through the aperture. 

It is also important to avoid wavefront distortions resulting from imperfect optical surfaces, such as 
reflective or transmissive substrates used in the fabrication of the LC-SLM. Analysis shows that 
imperfections of as little as L/20 can impose performance penalties of approximately 5 percent, and much 
more as the magnitude of surface imperfection grows. If those imperfections cannot be eliminated, 
wavefront correction is necessary, and the LC-SLM may be capable of performing those corrections. 

Further analysis was conducted for one limiting case of the LC-SLM, i.e., electrodes large enough 
that fringe fields and elastic behavior contribute very little to device performance. In that case, the 
phase retardation profile produced by an LC-SLM will appear to be stair-stepped. In order to maximize 
LC-SLM speed and minimize scattering losses, the device should be as thin as possible, which means 
removing the In degeneracies from the phase profile. That step produces phase “resets” in the final stair- 
stepped phase profile. 

It is already known that these resets do carry a penalty, but prior analyses considered only the case in 
which individual stair-step approximated blaze periods contain fixed and identical numbers of segments. 
This limits any segmented liquid crystal device to a finite number of discretely separated steering angles. 
That result has been extended to the more general case of phase resets falling on. and not just between, 
electrode segments. There is a correct choice for the phase value to be produced at electrodes for which a 
reset falls within bounds of that electrode. If the phase stair-step is obtained by averaging the value of the 
underlying blazed phase profile over each electrode before removing the 2n degeneracy, analysis shows 
that it does not matter where the reset falls. As a result, we can say that continuous beam steering is 
possible without any added penalty. 

As is already known, increasing the steering angle (i.e., blazes of smaller size) or increasing the 
number of waves being corrected by a given fixed aperture results in increased diffractive losses. Use 
of eight steps per blaze will result in a theoretical 5 percent loss, measured at the receiver's collecting 
aperture. 
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Incident light 

Figure 1 . — Approximate appearance of the director config- 
uration of three periods from a liquid crystal spatial phase 
modulator, with lower electrode voltages selected to 
produce blaze-phase grating. 


Integer 
multiples 
of 2-rr 



Figure 2. — Ideal linear phase ramp, characteristic of deflected 
plane waves, with its 2-u phase degeneracies removed. 


Phase retardation 
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(a) 


Position within grating 


Phase retardation 



Figure 3. — Phase retardation profiles produced by (a) large 
pixel stepping, and (b) small pixel OPA, showing evidence 
of “phase smoothing” resulting from increased significance 
of the spatial extent of fringing fields and liquid crystal 
elasticity. 
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Far field angle (.irad 

Figure 4. — Intensity profiles, in the far field, of a Gaussian beam which 
has been truncated in the plane of the minimum waist by an aperture. 
Intensity is normalized to that of an unhindered beam. 
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Figure 5. — Effect of sinusoidal surface undulations on far field inten- 
sity distribution of Gaussian beam reflected from a non-ideal mirror 
surface. Intensity is normalized to that of an unhindered beam. 
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Figure 6. — Strehl ratio versus surface undulation amplitude, for 3 
different undulation spatial frequencies, for a Gaussian beam 
reflected from a non-ideal mirror surface. 
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Figure 7. — Coverage of far field, from an OPA for which steering 
resolution is limited to the spot resolution. Not all parts of far 
field will be able to see the signal at its most intense. 
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Figure 8. — Ideal linear phase ramp, after application of 
modulo-2TT operation, overlaid onto patterned electrodes, 
for the case in which phase resets do not coincide with 
the boundary between electrodes. 



Figure 9. — For the case of Figure 7, applying the phase-averaging 
operation to the interperiod electrode, after first applying the 
modulo-2-rr operation to the ideal linear phase ramp. 
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Figure 10. — For the case of Figure 7, applying the modulo-2-n 
operation to the interperiod electrode, after first applying the 
phase-averaging operation to the ideal linear phase ramp. 
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Figure 1 1 . — Sections of grating phase profile, after (a) application 
of modulo-2u and then phase-averaging operation, and (b) 
reversing the order of those operations. 



Figure 1 2. — Strehl ratio versus number of electrodes 
per ideal blaze period, showing the effects of 
non-integer numbers for the two methods of 
assigning phase to interperiod electrodes. 
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Figure 13. — Idealized configuration of liquid crystal director, provid- 
ing a linearly increasing (l-r) optical path length for rays passing 
through the layer at normal incidence. This would constitute one 
period of a hypothetical liquid crystal blazed grating. 



Figure 1 4. — Gray scale map of the transverse field propagating through and away from 
a section of hypothetical liquid crystal grating constructed from the building block 
shown in figure 13.The light is normally incident on the grating from the bottom of 
the figure. Approximately 3 grating periods are shown in the lower 1/3 of the picture. 
The upper 2/3 shows the exiting light propagating in a vacuum. 
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Figure 1 5. — Similar to figure 1 4, but showing the phase component of the light. 



Figure 16. — Similar to figure 14, but showing the intensity of the light. 
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